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We discuss decoherence in charge qubits formed by multiple lateral quantum dots in the framework of the 
spin-boson model and the Born-Markov approximation. We consider the intrinsic decoherence caused by the 
coupling to bulk phonon modes. Two distinct quantum dot configurations are studied: (i) Three quantum dots 
in a ring geometry with one excess electron in total and (ii) arrays of quantum dots where the computational 
basis states form multipole charge configurations. For the three-dot qubit, we demonstrate the possibility of 
performing one- and two-qubit operations by solely tuning gate voltages. Compared to a previous proposal in- 
volving a linear three-dot spin qubit, the three-dot charge qubit allows for less overhead on two-qubit operations. 
For small interdot tunnel amplitudes, the three-dot qubits have Q factors much higher than those obtained for 
double-dot systems. The high-multipole dot configurations also show a substantial decrease in decoherence at 
low operation frequencies when compared to the double-dot qubit. 

PACS numbers: 73.21.La, 03.67.Lx, 73.23.Hk 



I. INTRODUCTION 

The realization of a solid-state qubit based on familiar and 
highly developed semiconductor technology would facilitate 
scaling to a many-qubit computer and make quantum com- 
putation more accessible. 1 The earliest proposal of a quantum 
dot qubit relied on the manipulation of the spin degree of free- 
dom of a single confined electron. 2 An attractive point of that 
proposal is the large spin decoherence time characteristic of 
semiconductors; a drawback is that it requires local control of 
intense magnetic fields. As an alternative, a spin-based logi- 
cal qubit involving a multiple quantum dot setup and voltage- 
controlled exchange interactions was devised, 3 but at the price 
of considerable overhead in additional operations. 

While spin qubits remain promising in the long term - note 
in particular several recent experimental advances 4,5 as well as 
further theoretical development of multiple-quantum dot spin 
qubits 6 - charge-based qubits in quantum dots, in analogy to 
superconducting Cooper-pair box devices /- 8,9 ' 10 are also wor- 
thy of investigation. Employing the charge degree of free- 
dom of electrons rather than their spin brings a few impor- 
tant practical advantages: No local control of magnetic fields 
is required and all operations can be carried out by manip- 
ulating gate voltages. The simplest realization of a charge 
qubit is a double quantum dot system with an odd number of 
electrons / W 3 ' 14,15,16 ' 17 One can view this system as a dou- 
ble well potential: The unpaired electron moves between the 
two wells (i.e., quantum dots) by tunneling through the po- 
tential barrier. The logical states |0), |1) correspond to the 
electron being on the left or right. The barrier height deter- 
mines the tunneling rate between the dots and can be adjusted 
by a gate voltage. The resulting bonding and antibonding 
states can also be used as the computational basis. Recently, 
three groups have implemented the double-dot charge qubit 
experimentally ! 18 i 19 i 20 

Charge qubits are susceptible to various decoherence mech- 
anisms related to charge motion. Strong damping of coher- 



ent oscillations was observed in all quantum dot experimental 
setups / 8 ' 19 ! 20 with quality factors in the range 3-10. Note that 
a change in the state of the qubit involves electron motion be- 
tween quantum dots, which can in general couple very effec- 
tively to external degrees of freedom such as phonons, charge 
traps in the substrate, and electromagnetic environmental fluc- 
tuations. These noise sources lead to decoherence times much 
shorter than those observed in spin qubit systems. Thus, one 
is tempted to try to find new setups where oscillations between 
qubit states involve a minimum amount of charge motion. For 
instance, in qubits based on multiple quantum dots one can 
pick logical states where charge is homogeneously distributed 
in space. Another approach is to create a multiple dot structure 
with symmetries that forbid coupling to certain environmental 
modes within the logical subspace. 21 A recent attempt along 
this direction is found in Ref.l22l 

In this paper we argue that it is not generally possible to 
avoid decoherence in multiple-quantum-dot charge qubits by 
simple geometrical constructions. The spreading of charge 
uniformly over a multiple-quantum-dot logical qubit does not 
avoid decoherence. However, the coupling to bosonic envi- 
ronmental modes, such as phonons and photons, can be very 
substantially attenuated in some circumstances. 

In order to demonstrate these assertions, we analyze in de- 
tail two prototypical extensions of the double-quantum dot 
charge qubit. We first consider a qubit consisting of three 
quantum dots forming a ringlike structure and only one ex- 
tra electron, as shown in Fig.Q] Multiple-quantum-dot qubits 
with a ringlike structure resemble a proposal by Kulik et alP 
to use persistent current states in metallic rings for quantum 
computation. Unlike the double-dot qubit case, the ground 
state in a three-dot qubit can be truly degenerate with corre- 
sponding wave functions having a uniform charge distribu- 
tion. At first, this raises the hope that decoherence mech- 
anisms involving charge inhomogeneities (such as phonons 
or charge traps) would be inhibited due to mutual cancella- 
tions. However, we shall see below that the computational 
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FIG. 1: (Color online) Schematic illustration of a three-quantum-dot 
qubit with only one extra, unpaired electron. The external tuning 
parameters are the strength of the tunneling couplings (vi, V2, and 
V3) and the magnetic flux $ = <j>\ + <j>2 + <f>3 through the qubit. The 
latter is used solely to define the working point of the qubit. 
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FIG. 2: (Color online) Eigenenergies of the three-dot qubit as func- 
tion of the magnetic flux. The working point at $/3 = <f) — w P er 
bond is indicated by the arrow. At this point, clockwise and coun- 
terclockwise persistent current states are degenerate, and the charge 
distribution is homogeneous throughout the space spanned by the 
computational basis. 



basis states can be distinguished by phonon and electromag- 
netic baths through the electron phase variations along the 
ring. That, in turn, leads to dephasing and decoherence. This 
problem is intrinsic to all quantum-dot-based charge qubits. 
Nevertheless, the Q factor in these three dot qubits can be 1-2 
orders of magnitude larger than in the corresponding double- 
dot qubits, a substantial improvement in coherence. 

Second, we show that planar quantum dot arrays in the form 
of high-order multipoles can be more efficient in reducing the 
coupling to acoustic phonons in multiple-quantum dot qubits. 
Our work extends and analyzes in detail a recent proposal to 
create a decoherence-free subspace with charge qubits. 22 

While it is well known that condensed-matter environments 
tend to produce time and spatial correlations in their interac- 
tion with qubits, 24 here we assume that the Markov approx- 
imation provides reasonable estimates of decoherence rates. 
In particular, we employ the Redfield equations in the weak- 
coupling, Born-Markov approximation to describe the time 
evolution of the reduced density matrix of the qubit system. 25 

The paper is organized as follows. In Sec.[n]we study in de- 
tail a three-dot charge qubit. Single- and two-qubit operations 
are presented, as well as the coupling to a bosonic bath. We 
consider in detail the particular case of acoustic piezoelectric 
phonons, which is relevant to III-V semiconductor materials 
at low temperatures. Also in Sec. [II] we evaluate decoherence 
and energy relaxation rates using the Redfield equation for- 
malism. In Sec. [HI] we present a multiple-quantum dot logical 
qubit structure that minimizes the coupling to environmental 
modes which couple to charge. We also analyze phonon de- 
coherence in these systems and compare Q factors with those 
obtained with double-dot charge qubits. Finally, our conclu- 
sions are presented in Sec.[lV] The Appendix contains math- 
ematical details of the two-qubit operation of the three-dot 
charge qubit of Sec. [TTJ Throughout this paper we assume 
H = 1 and ks = 1 ■ 



II. THE THREE-DOT CHARGE QUBIT 

A simple example of a multiple-dot qubit with charge der- 
ealization consists of three quantum dots in a ringlike geome- 
try, as shown in Fig. \T\ In practice, this system is created by 
laterally confining electrons in a two-dimensional plane; the 
confinement is electrostatic, controlled through electrodes sit- 
ting above the plane. Consider gate voltages on the electrodes 
such that the three dots share one excess, unpaired electron, 
while all configurations with a different number of excess 
electrons become energetically inaccessible due to the large 
charging energy of the doXsM The spin degree of freedom is 
not relevant for our discussion and electrons will be assumed 
spinless unless otherwise specified. Thus, the system lives in 
a three-dimensional Hilbert space. The electron can hop be- 
tween dots through tunneling. The tunneling matrix elements 
and the on-site energies are controlled by the gate voltages. As 
will be clear shortly, it is convenient to apply a weak magnetic 
field perpendicular to the plane containing the dots. 

The three natural basis states place the electron on dot A, 
B, or C: 

\A) = c A |vac), |S)=4|vac), \C) = c c |vac), (1) 

where cj, are creation operators and |vac.) is a reference state 
where all dots have an even number of electrons. In this basis, 
the Hamiltonian takes the matrix form 





( E A 


-vi e^ 1 


-v 3 e-^ 3 




-vi e-** 1 


E B 


-v 2 e^ 2 




\ -u 3 e^ 3 


-v 2 e-^ 2 


E c 



where Ea,Eb, and Ec are the on-site energies, Vi are the 
tunneling strengths between pairs of quantum dots, and 4>\ + 
02 + ^3 = $ is the total magnetic flux through the ring. 
Let us specify the qubit by setting vi = v 2 = V3 = v > 0, 
Ea — Eb = Ec = 0, and 4>i — 4> 2 — 4>z = $/3 = tt. In this 
configuration, two degenerate eigenstates |+) and |— ) have 
the lowest energy, E± = —v (Fig. [2]). They carry clockwise 
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and counterclockwise persistent currents and form the compu- 
tational basis. 27 The third, excited, eigenstate |T) has energy 
E e — 2v and is current-free. The eigenvectors are 

\T) = -±=(\A) + \B) + \C)), (3) 

1+) = -^(\A)+e^\B) + e-^\C)), (4) 

h) = -^(\A)+e-^\B)+e^\C)), (5) 

with (3 — 2ir/3. Clearly, the charge distribution is spatially 
uniform for all three states. 

It is worth noting that the topology of the three-dot qubit 
and its use of persistent currents of opposite direction as log- 
ical states closely resemble the Josephson persistent current 

I 

/ 2v+l(5 1 + 5 2 + 5 3 ) 
H = -i (fce* + S 2 + 5 3 e- i/3 ) 
\-ils ie -i0 + S 2 + S 3 e i P) 



The computational subspace corresponds to the lower-right 
2x2 block. Evidently, we stay within the computational sub- 
space as long as 5\ = 5-2 = S3. However, this also implies that 
there is no coupling between the computational basis states 
|+) and |— ). For 8ie % @ + 62 + ^e" 1 ^ 7^ 0, coupling within 
the computational subspace is possible, but there is a finite 
probability of leaking out into the state \T). The leakage can 
be kept small as long as v 3> |#i,2,3|. Alternatively, one can 
incorporate the third level into the single-qubit operations, as 
in Ref. 23. For the following case study, we assume that the 
leakage from the computational subspace is negligible. 

Using the Pauli matrices ax, cr 2 , and 03, as well as the iden- 
tity matrix a n , we can express the Hamiltonian in the compu- 
tational basis in terms of a pseudospin in a pseudomagnetic 
field h plus a constant, 

H s — E <7 + h x (7i + h v (7 2 + h z cr 3 , (7) 

where E$ = —v — (Si + 82 + S3) /3 and 



2 fx ft + M 


(8) 


3 V 2 2 )■> 


Si - S3 


(9) 


V3 ' 


— Vif/VS. 


(10) 



We only need to vary two out of the three pseudomagnetic 
field components in order to perform single-qubit operations. 
Thus, we can operate the qubit at constant magnetic flux (and 
set ip = 0, h z = 0) and vary only the Si via gate voltages. 
If we furthermore fix the coupling V2 = v, 62 = 0, we find 
that the qubit is controlled by the sum and difference of the 



qubit studied in Ref. 28 or the proposed atomic Josephson 
junction arrays. 29 However, the similarities stop here as the 
underlying physics is very different. We will focus our dis- 
cussion on the quantum dot charge qubit case only. 



A. Single-qubit operations 

In order to be able to perform quantum gate operations, we 
have to allow for deviations from the degeneracy point. This 
is done by varying the tunneling coupling and/or the magnetic 
flux. It is convenient to introduce the (small) parameters Si, 
S2, S3, and ip such that vi = v + Si, v 2 = v + 82, V3 = v + S3, 
and (p <C 1 with <p — $ — 3it. To linear order and using a 
{|T), |+), |— )} basis, we find that the Hamiltonian expanded 
around the degeneracy point can be written as 



-l{Sie^ + S2 + S 3 e-^) \ 
§ (Sie-V + 8 2 + 8 3 e^) . (6) 

-v + vip/V3 - § (Si + 8 2 + S 3 ) J 

I 

variation of two intra-qubit couplings, h x oc (Si + S3) and 
h y oc (81 — S3), that can be adjusted by tuning the respective 
gate voltages around the symmetry point. 

B. Two-qubit operations 

In order to perform two-qubit operations, such as the SWAP 
or CNOT gate, we have to couple two three-dot qubits (called 
I and II hereafter). In principle, this can be done in either a 
tip-to-tip or base-to-base coupling scheme, as shown in Fig. [3] 
Since the number of excess electrons in the composite system 
is equal to two, states where two electrons occupy the same 
qubit have to be included in the basis of the two-qubit Hilbert 




FIG. 3: (Color online) Possible implementations of two-qubit gates 
using three-dot qubits. (a) Coupling via a single dot (tip-tip geome- 
try); (b) coupling via two dots (base-base geometry), (c) A possible 
implementation of a qubit chain in the base-base configuration. 



-§ (£ie-^ + £ 2 + <5 3 e i/3 ) 
-v-vtp/V3- I (61+ S2+S3) 
l(Sie^ + 5 2 + 8 3 e-^) 
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space. The basis of the two-qubit Hilbert space reads thus 



11) 


= l+>il+>ii, 


(ID 


|2> 


= |+>i|-)n, 


(12) 


|3) 


= |-)i|+)n, 


(13) 


1 -1 \ 

I 4 ) 


1 \ 1 \ 


(14) 


I 5 ; 


= c ai c bi I vac.), 


(o) 


1 6) 


t t i \ 


(16) 


1 7) 


= c m c ci l vac ->> 


(17) 


IS) 


= C AII C SII vac -)> 


(18) 


|9) 


= c ylII c CII vaC -)' 


(19) 


1 10) 


= C BII C CII vac -)- 


(20) 



Here, two types of states have been neglected: First, states 
with double occupancy of a single dot since the charging en- 
ergy is assumed to be very large. Second, although the \T)\ 
and \T)u states couple to the double-occupied states |5) to 
1 10) through the inter-qubit hopping terms, they are gapped by 
an energy of order v, which is assumed much larger than the 
effective two-qubit interaction amplitude t' 2 /Ui (see below). 
Therefore, they were not included in the two-qubit Hilbert 
subspace^ 

The Hamiltonian for the inter-qubit interaction in the tip-tip 
setting shown in Fig. [2a) reads 



tip 

I— II — 



-t'{c\ 



BI °CII 



-en 



c Bl)- 



(21) 



Similarly, the base-base coupling presented in Fig. [2b) is gov- 
erned by the Hamiltonian (see also the Appendix) 



rrbasc 
^1 II 



-t'(c 



Bl C CII 



J CII ^Bl) 



t"(c 



BII 



-fill CI/ 



(22) 

where we have chosen the gauge for the vector potential as- 
sociate to the perperdincular magnetic field to be parallel to 
the inter-qubit tunneling paths. We assume that the inter-qubit 
tunneling amplitudes t' and t" satisfy < t' , t" <C v <C Ui, 
where Ui is the inter-dot charging or capacitive coupling en- 
ergy (i.e., the change in the energy of one dot when an electron 
is added to one of the neighboring dots). In other words, the 
capacitive coupling between dots must be sufficiently strong 
so that states with two or zero excess electrons in a qubit are 
forbidden. Due to the proximity between dots of neighboring 
qubits, some small inter-qubit capacitive coupling will also 
exist. Although we will neglect such coupling in the discus- 
sion below, these additional charging energies can be included 
without substantially modifying our results. In particular, we 
note that the inter-qubit capacitive coupling does not interfere 
with single-qubit operations. Note also that the presence of 
a magnetic flux requires the dots A, B, and C to be always 
arranged in a clockwise order. 

Next, the large charging energy separation between the 
single-occupancy states |1) to |4) and the double-occupancy 
states |5) to 1 10) allows us to separate the two-qubit compu- 
tational subspace from the rest of the Hilbert space. In order 
to do so, we use a Schrieffer- Wolff transformation, 31 which 



amounts to a second-order perturbative expansion of the ef- 
fective Hamiltonian in the ratio of the inter-qubit tunneling 
magnitude to the charging energy. To this end we insert the ex- 
pressions for | +) and | — ) from Eqs. ©-(O into Eqs. (fm>-(fT4l> 
and express the computational basis states |1) to |4) in terms 
of creation operators acting on the vacuum state. Further, us- 
ing the basis vectors in Eqs. (fTTT > to ( 1201 ), one can easily com- 
pute the full six-dot Hamiltonian in the basis of states |1) to 
1 10). Noting that one can obtain the tip-tip Hamiltonian from 
the expression for the base-base case by setting t" = 0, we 
evaluate the more general case of the base-base coupling, see 
Fig. [23 and Eq. t22i . The details of the computation, i.e. the 
full matrix representation of this Hamiltonian, as well as its 
reduction to the two-qubit computational basis by performing 
the Schrieffer- Wolff transformation, are shown in Appendix 
lAl The result for the reduced Hamiltonian takes a rather com- 
pact form which, for the tip-tip case, reads 



frtip _ 

^i-ii — 



( 



\ 



4 

-i/3 

-2 



4 

-2e*P 

-if3 



e i0 

4 



-2 

e -i0 
4 



(23) 



Note that this reduced Hamiltonian acts on the subspace 
formed by the states {|1), . . . |4)} defined in Eqs. (Qj} to lfl4l . 
Up to the common prefactor —t' 2 / (9 Ui), the eigenvalues of 
H^ u are E 1 = 0, E 2 = 4, E 3 = 6, and E 4 = 6, with the 
respective eigenvectors equal to 



\Ex) 
\E 2 ) 
\Es) 



i(|l)- e ^2)- e -^|3> + |4)), 

i(|l> + e tf|2) + e-*|3> + |4) 
1 



-iff 



|3) 



(24) 
(25) 
(26) 
(27) 



The critical question now is whether this setup permits a 
convenient two-qubit operation, such as a full SWAP. It is 
straightforward to show that the answer is positive, even in 
the simple tip-tip coupling scheme. To see that, suppose we 
initialize the qubits in state |2) and now search for the time t 
after which the qubits have evolved onto the (swapped) state 
1 3) under the action of Hj^ u . The square of the resulting con- 



dition, 



(3|e- 



|2) 



= 1, is readily evaluated and yields 

t s = tt/2 [t' 2 /{9 Uj)]^ 1 as the (shortest) time for which the 
tip-tip coupling t' has to be turned on in order to implement 
the SWAP gate. 

For a comparison with the (linear) three-dot spin qubit 
scheme proposed by DiVincenzo et a/., 3 let us briefly discuss 
the implementation of the CNOT quantum gate. A CNOT 
can be done straightforwardly using two VSWAP operations 
(SWAP gates of duration ts/2) and seven one-qubit gatesr^ 
e.g., by utilizing the scheme in Ref.[32j. Consequently, we 
find that the realization of one- and two-qubit operations for 
the present three-dot charge qubit is considerably simpler than 
for the proposal by DiVincenzo et al. where many more steps 
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were necessary to implement a CNOT. One reason is the com- 
plexity of the one-qubit rotations - for the logical spin-qubit, 
one-qubit operations alone require three spin exchange inter- 
action pulses. For the CNOT gate, this implies at least 19 
pulses with 1 1 different operation times. Compared to the 9 
pulses needed for the three-dot charge qubit, the practical ad- 
vantages of the qubit and computation scheme proposed here 
are evident. 



C. Coupling to a bosonic bath 

The charge qubit couples to a variety of environmental de- 
grees of freedom. We study in particular the decoherence 
caused by gapless bosonic modes that sense charge fluctua- 
tions in the dots, such as phonons. We assume that all quan- 
tum dots couple to the same bath. The Hamiltonian describing 
the non-interacting bosonic modes in this case is 

H B = Y J ^b\b^ (28) 
q 

with q denoting the boson linear momentum and aj q its disper- 
sion relation. The coupling between the dots and the bosons 



is assumed to be governed by the bilinear Hamiltonian 3 

oson — 

(29) 

where Nk is the number operator of the /cth dot, and 

« t =A q fe ,R 'r (30) 

Here, A q represents the electron-boson coupling constant and 

Pq and Rfe are form factor and position vector of the kth 
dot, respectively. Note that all geometrical information is 
contained in the coefficients a^. Since we have exactly 
one excess electron on the three-dot system, the constraint 
Na + Nb + Nc = 1 must be satisfied. Therefore, the system- 
bath Hamiltonian in the basis {\A), \B), |C)} reads 

( OLA \ 
#sb = V a B (6 q + 6_q). (31) 
q V a c J 

Projection of this Hamiltonian onto the subspace spanned by 
|+) and |— ) defined in Eqs. (01 and <(5j constrains the coupling 
to that subspace, yielding 



q L 



ola ^ ) CTl 2~ _ a 172 



(32) 



where a term proportional to <jq has been dropped. The pres- 
ence of two terms with different functional dependence on q 
indicates the coupling to two bath modes, which will be de- 
noted by the indices 1 and 2 in the following. There would be 
a third bath mode, proportional to er 3 , if the charge distribution 
were not the same for the two logical states. The advantage 
of having a homogeneous charge distribution for both states 
in the computational basis, leading directly to the cancella- 
tion of this third mode of decoherence, is evident here. It is 
important to remark that charge homogeneity can be achieved 
without the assumptions of homogeneous tunneling or equal 
capacitances: as long as one can tune the gate voltages in the 
quantum dots independently, one can arrange to have one ex- 
tra electron equally shared among the three dots. 

It is convenient to rewrite the system-bath Hamiltonian in 
the standard spin-boson formal 



describe the system part and the corresponding bath part is 
given by $ 1>2 = J2 q 9^ (& q + &-q). with 



9? 

,(*) 



2a a — as — etc, 
a B - ac ■ 



(35) 
(36) 



Assuming all P q to be the same, the following relations 
among the a^ can be obtained: 



OLA = A q Pq, 

a B 



a c = a A e l(Rc -' RA) - cl = aAe lT > c 



(37) 
(38) 
(39) 



where the last two equations define the phases t\b and rjc- 
This completes the specification of the qubit-bath coupling. 



fr Ty~ tf\ i #f» D. The Redfleld equation 

-"SB = + A 2 V 2 , (JJ) 

We now investigate the qubit decoherence due to the 
where bosonic bath by determining the time relaxation of the sys- 

tem's reduced density matrix. We use the Born and Markov 
_ approximations and the Redfield equation. 25 In this formalism 

Ki=a\/& and Ki = — <72/2\/3 (34) the reduced density matrix of the system (qubit) is obtained by 
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integrating out the bath degrees of freedom and assuming that: 
(i) the coupling to the bath is weak, so leading order pertur- 
bation theory is applicable (the Born approximation), and (ii) 
the bath correlation time is much shorter than the typical time 
scale of operation of the qubit, so that system-bath interaction 
events are uncorrelated in time (the Markov approximation). 

The time evolution of the reduced density matrix is given 
by the Redfield equation?5i2£ 

p(t) = -i[H s (t),p(t)] (40) 

+ J2 { [K(t)p(t),K a ] + [K a ,Ai(t)p{t)] }, 
0=1,2 

where the time-dependent auxiliary matrices A a (t) which en- 
code the bath correlation properties are defined by 



,9=1,2 



Hs(t) 



(41) 



The thermal-average bath correlation functions, 

B a p(t) = * o (0)>, (42) 

can be written in terms of spectral functions, 

^M^E^^^-^q), (43) 



and the dispersion relation w q are both isotropic. Second, let 
us assume that the electronic density in the dots has a Gaussian 
profile, p(r) = S(z) e~ r ^ 2a ^/(2na 2 ), resulting in 



J d 3 rp{r) e~ iqr = e -( Q « sin9 ) 2 / 2 , (48) 

where (q, 9, ip) are the spherical coordinates of the boson 
wave vector. Then, the three-fold symmetry in the plane 
causes V\ 2 (ui) to vanish and 



(49) 



with 



fig 2 
2^2 



du>r. 



dq 



tt/2 



dti 



x [1 - J {qD sin 6)] 



(50) 



where u> = u> q , ft is the crystal unit cell volume, and D is the 
distance between dots. 

For III-V semiconductor materials at low temperatures, 
the most relevant bosonic modes are piezoelectric acoustic 
phonons, 37 for which we have X q — xswg^hJqTl and uo q — 
sq. Here, <? p h is the dimensionless electron-phonon coupling 
constant and s is the phonon velocity (for GaAs, g p h ~ 0.05 
and s w 5 x 10 3 m/s)^ 



and the boson occupation number n B (oS) = (e"/ T — 1) 1 : 



B af3 {t) = / dwv af} {uj) (44) 
Jo 

x{e wt n B {uj) + e-^ t [\ + n B {u)])} . 
Performing the sum over q in Eq. ( |43l , we find 

V n = 2^ |A q Fq| 2 5{ljJ - UJg) 

q 

x [3 - 2(cos7?b - cos -qc) + cos(?7 B - Vc)] (45) 

^22 = 2^|A q P q | 2 5(w-w q )[l-cos(7/ B -7 ?c )] J (46) 
q 

V\2 = 2^ |A q F q | 2 6(U - UJg) 



-"IB _ p—lTIG 



+ ism(r) B - r)c)] , (47) 



with j/21 = ^12- When the bath is sufficiently large, the sums 
over the vector q in Eqs. (|43T>-(|47]| can be converted into three- 
dimensional integrals. 

A few simplifying but realistic assumptions can be made at 
this point. Let us first assume that the coupling constant A q 



E. Decoherence rates 

We now solve the equation-of-motion for the reduced den- 
sity matrix explicitly for a case in which the decoherence rate 
can be obtained directly. Consider a constant pulse applied to 
the qubit at t = such that h y = h z = Q and h x = A > 0. For 
t > 0, the A a matrices are constant and given by 



Ai = 7oCi/2 , 

A 2 = -(1/2V3) (7^2+7^3) 
The (complex) relaxation rates are given by 



70 ee / dtB 22 (t), 
10 



dtB 22 {t) sin(2Ai). 



The relaxation part of Eq. (l40l then reads 



(51) 
(52) 



(53) 



7c = / dtB 22 (t) cos(2At), (54) 
la 



(55) 
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2^ IL^p.Aj+n.c.) 6 ^ 12 -^ 2 pn -P22J 6 i -P12 - P12 Pn-P2 2 J 



) 




'0 


0) 


- P12 y 
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(56) 



where the single and double primes denote real and imaginary 
parts, respectively. They can be easily evaluated, yielding 



7o 



0. 



^ = -i/(2A)coth - , 



A 



T 



(57) 
(58) 



7 S = 



^*,(2Ay)coth(^y (59) 



P12 ~ P12 P22 - P11 

Pll - P22 Pl2 - P\ 2 



The Liouville term in Eq. d40b is obtained straightforwardly 

-i[Hs,p] = -iA[a 1 ,p\=-iA 

Introducing Eqs. (T56b and ( f60b into d40b . we obtain 
P11 = -2 I A 



(60) 



7s 
G 



• / 7c 1 1 7c 

Pl2 — Pl2 T "g"j 



P12 + y (1 - 2pn) 



P12 



-A(l-2p u ) 



(61) 

(62) 
(63) 



where we have split the off-diagonal term pi 2 into real and 
imaginary parts, p[ 2 + «pi' 2 - 

In order to identify energy and phase relaxation rates, we 
rewrite the elements of the reduced matrix in the eigenbasis 
of the system Hamiltonian, 



\E = ±A) 



1 

71 



l+)±l-». 



resulting in 



P11 



P12 



P12 



T Pll + y 



- ^2A 
2bp{ 2 . 



P12 



3P12: 



(64) 

(65) 

(66) 
(67) 



The solution of the diagonal term is straightforward, 

Pu(*) = l + [#i(0)-l] e~^ t/3 , (68) 
which allows us to read directly the energy relaxation time, 

(69) 



7c 



For the off-diagonal term, one finds that the real part is given 
by 



P12W =Pi' 2 (0)e- t/T2 cos(w c t), 



(70) 






100 125 
v [ueV] 



200 



FIG. 4: (Color online) Comparison between the Q factors of a three- 
dot and a double-dot charge qubit coupled to piezoelectric acoustic 
phonons. The parameters used are: a = 60 nm, D — 180 nm, 
s = 5 x 10 3 m/s, T — 15 mK, and 3 P h = 0.05, which correspond 
to realistic lateral quantum dot systems in GaAs. Here the variable 
v denotes the interdot tunnel amplitude. Note that for double dot 
qubits, A = v, while for three-dot qubits we assumed A = 0.1 v. In 
the inset we show the same Q factors when the oscillation frequency 
(rather than v) is fixed. In this case the curves only differ by a factor 
of3. 



where the phase relaxation time is equal to 



To 



G 



and the frequency of quantum oscillations is given by 



UJr 



'2A 2A 



7^ 
36 



(71) 



(72) 



Note that T 2 = 22\, as well-known for the super-ohmic spin- 
boson model in the weak-coupling regime322 

Except for a factor of three in the relaxation rates, Eqs. 
(|69]l-(|72l are identical to those found in Ref.[l5| for a double- 
dot charge qubit. However, one has to recall that while in the 
double-dot qubit A is the interdot hopping matrix element v, 
for the three-dot qubit it takes a much smaller value, of the 
order of <5i 2,3- The decoherence times will be longer for the 
three-dot qubit, but so will be the quantum oscillation period 
and the single-qubit gate pulses. Therefore, it is meaning- 
ful to compare the quality factor of the the three-dot qubit to 
that obtained for the double-dot qubit for a fixed magnitude 
of v, which is a common experimental parameter to both se- 
tups. The comparison for the case of piezoelectric acoustic 
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phonons and realistic GaAs quantum dot geometries (data for 
the double-dot qubit was obtained from Ref. 15) is shown in 
Fig- E] The Q factor is defined as 



We assume ui c w 2 A since A 3> j' c , j' s in the weak-coupling 
regime. 

The improvement in the Q factor is substantial for small 
tunnel amplitudes. A similar result was previously found by 
Storcz et al. in Ref.|40j when considering the phonon-induced 
decoherence in a system of two double-dot charge qubits with 
a small tunnel splitting ("slow tunneling"). There, the dom- 
inant quadrupolar contribution to the two-qubit decoherence 
yields a uj^ dependence for the Q factor. In our case, the extra 
protection in the three-dot qubit compared to the slow tunnel- 
ing double-dot system arises mainly because the oscillation 
frequency cu c (i.e., the amplitude of the transverse pseudomag- 
netic field) is smaller in the three-dot qubit by the ratio A/v 
[see Eq. ©]. This ratio must be kept small in order to avoid 
leakage from the computational basis. In Fig. [4] it was set to 
0.1. However, for a fixed oscillation frequency (see inset in 
Fig. @), the Q factors for these two qubits differ by only an 
overall factor of three. 

To summarize up to this point, our study indicates that using 
a computational basis with a homogeneous charge distribution 
improves the quality of the qubit but does not rid it from de- 
coherence completely. The reason lies in the fact that bosonic 
modes propagating in the xy plane can pick up distinct phase 
shifts when interacting with different dots [see Eqs. ( f3Tb to 
d39l>l. However, there is no complete destructive interference 
along any direction of propagation in the plane, as can be seen 
from Eqs. ( f33T > and ( f36b . In fact, one can show that the same 
is true for any ringlike array of dots that share a single excess 
electron. 



III. CHARGE QUBITS IN MULTIPOLE 
CONFIGURATIONS 

As recently proposed by Oi et al., 22 there is another way 
in which the geometry of the quantum dot qubit array and its 
charge distribution can be chosen to minimize the coupling to 
environmental degrees of freedom. Here we demonstrate how 
their idea can be extended to multiple-dot charge qubits cou- 
pled to gapless bosonic modes. It turns out that by reducing 
the computational space to particular multipole charge con- 
figurations one can substantially reduce the coupling to bath 
modes at low frequencies. We consider qubits and basis states 
as shown in Fig. [5] The qubit consists of a planar array of dots 
with alternating excess charge. Note that the operation of such 
a qubit is straightforward: The excess charge is only allowed 
to hop between every other pair of neighboring dots, namely, 
between dots numbered 2n — 1 and 2n, with n = 1, . . . , 2 P_1 , 
where p is the multipole order, I = 2 P (see Fig. 0. Tunnel 
barriers between alternating pairs of dots must be maintained 
small and fixed (to avoid leakage), while the remaining bar- 
riers have to be modulated in time to implement an X gate. 



1=2 1=4 1=8 




FIG. 5: The three lowest multipole charge qubit configurations 
(dipole, quadrupole, and octopole). The two computational basis 
states, jO) and |1), are indicated for for each configuration. Empty 
(filled) circles correspond to empty (occupied) quantum dots. The 
arrows indicate the pairs of quantum dots where excess charge can 
hop. 

The Z gate is implemented by inducing a small bias between 
even- and odd-numbered dots. Two-qubit operations can be 
implemented in analogy to the procedure discussed in Ref.l36l 
in the context of Ising-interaction based two-qubit operations. 

The basis states for each multipole configuration have com- 
plementary charge distributions that tend to cancel out the 
coupling to phonon modes propagating along certain direc- 
tions in the xy plane. The number of such directions in- 
creases with the multipole order, resulting in an attenuation 
of the overall coupling to phonons at low frequencies (large 
wavelengths). The crossover frequency where this attenua- 
tion occurs is Across ~ s/di, where di is the radius of the 
dot array. At high frequencies, however, when the phonon 
wavelength is much smaller than the radius di, decoherence 
becomes stronger because phonons can resolve the internal 
structure of the qubit and disturb charge motion between indi- 
vidual pairs of dots. 

In order to demonstrate these effects, let us derive an ex- 
pression for the spectral function of the qubit-bath system. For 
simplicity, we assume that all dots in the qubit are identical. 
In this case, the bath modes couple to charge variations in the 
dots according to the Hamiltonian 

H SB = J2I2 akNk ( b l + b -*)' (74) 

q fc=l 

where I = 2 P and N). is the excess charge in the fcth dot. For 
the case of acoustic phonons, the coefficients ct^ were defined 
in Eq. ( 13 OI . Projecting this Hamiltonian onto the computa- 
tional basis (as shown in Fig. [5j, we find that, up to a constant 
term, 

H SB = K<P h (75) 
where K = —a z /2 acts on the qubit space and 

$1 = ^4')^ + ^), (76) 
q 
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acts on the phonon bath, with 



7« = £(-l) fe "* 

k=l 



I 

fe=l 



(77) 



It is convenient to choose the position vectors of the dots as 
R-fc = di(x cost^fc + y sinifik), where ipt- = (2iv/l)(k — 1) 
and di is the array radius: di = D/2 sm(iv/l), where D is the 
distance between neighboring dots. This yields 



iA q ^i 2 E (-v™ 



(78) 



k,j=l 



x exp 



2idiq sin sin (^p — 



sm 



<Pk - V] 



where (q, 9, ip) are the spherical coordinates of the wave vec- 
tor q. It is not difficult to see that g q = for 9 = ir/2 and 
ip = (2m — with m = 1, . . . , I. 

The spectral function can now be obtained in analogy to the 
calculation shown in Sec. Ill DI For a thermal bath of acoustic 
piezoelectric phonons, we find 



E K° 

q 




S(UJ- UJq) 
tt/2 / 

d9 sin# exp I — 
(-D'/Vo^sin, 



(79) 



+2 e (-ir Jb 

m— 1 



2<i;W . 

sin W sin . 

s V Z 



/m7r 



Implicit in Eq. ( f79l are the assumptions of in-plane isotropy of 
A q , -Pq, and o> q . Note that for I = 2 one recovers the spectral 
function for a double dot qubit obtained in Ref.|l_5|. The low- 
frequency behavior of the spectral density becomes apparent 
when we expand the Bessel functions in a power series, re- 
sulting in 



7r/2 



2 jo 

y (-i) fc 



d0 sin exp 



sin 2 6> 



sm I 



2 A- 



(80) 



where 



,(0 _ 



1/2— 1 

(-1)'/ 2 + 2 E (-l) m sin 2fc 

m— 1 

,(0 - 



(81) 



It is possible to show that a k = for k < 1/2 when / is an 



integer power of 2. Therefore, — > 0) 



,1+1 



For large i, 



this amounts to the appearance of a pseudo gap in the spectral 
function at low frequencies. The asymptotic behavior of the 
spectral function at high frequencies is also straightforward to 



obtain: One finds vi(uo — > oo) ~ l/u>. Thus, the tail of the 
spectral function raises with increasing multipole order. 

The structure of the computational basis is simple enough to 
allow for the qubit to couple to just one bath mode (in contrast 
to the three-dot qubit, where two modes couple to the qubit). 
Thus, the standard expressions for the relaxation times in the 
spin-boson model can be used. 39 The result is 



v(0 



-vi(2v)coth{- 



(82) 



where v is the inter-dot tunnel amplitude. Note that Eq. 
reduces to the result found in Ref.[l5| when / = 2, as expected. 
Provided that v is sufficiently smaller than the temperature, 
>yW v l _ Thus, by increasing the order of the multipole and 
maintaining a low frequency of operation, one can decrease 
the qubit relaxation rate by orders of magnitude without af- 
fecting the frequency of quantum oscillations. In Fig. [6] we 
show the Q factor of several multipole charge qubits as a func- 
tion of the inter-dot coupling v. Note that at low frequencies 
high quantum oscillations are much less damped for high mul- 
tipole configurations. This translates into single-qubit gates of 
much higher fidelity. Clearly, this gain in the Q factor has to 
be contrasted with the high complexity of operating a logical 
qubit comprised by a large number of quantum dots, as well 
with the slowness in operation. As the gating becomes slower, 
other sources of decoherence may become more relevant. 

It is also important to remark that, in practice, the pseudo- 
gap width, across, will shrink with increasing multipole order. 
This is because the dot array radius scales as di k,ID/2it for 
I 1. Therefore, for a fixed value of D, one has across ~ 
2irs/(lD) for I 3> 1. Finally, we note that the results dis- 
cussed above are valid for any gapless bosonic bath. Different 
dependences on q for the coupling constant \ q and disper- 
sion relation uj q will only change the power of the frequency- 
dependent prefactor in Eq. (f80b . 



— dipole 

--■ quadrupole 

- octopole 
hexadecapole 




FIG. 6: (Color online) Q factors for multipole charge qubits 
(I = 2, 4, 8, 16) coupled to piezoelectric acoustic phonons: Qi = 
lJc/kj^, where uj c » 2v [for 7"% see Eq. d82H. Physical and 
geometrical parameters are the same as those used in Fig. [4] In par- 
ticular, note that the inter-dot distance is fixed, D — 120 nm, for all 
configurations. 
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IV. CONCLUSIONS 

In this paper we have shown that, whereas there are no 
simple ways to completely protect charge qubits based on 
quantum dots from decoherence by gapless bosonic modes, 
a homogeneous charge distribution throughout the qubit is the 
most advantageous setup and provides the best possible pro- 
tection against decoherence. This result applies not only to 
the charge qubits in semiconductor-heterostructures that we 
focused on here, but, in principle, to charge qubits in gen- 
eral. Whereas certain aspects of the discussion need to be 
changed for, say, self-assembled quantum dots, single-donor 
charge qubits, 41 or Si-based quantum dot structures,— this 
does not affect the universal mechanism underlying our cen- 
tral result, namely that a specific (homogeneous) charge dis- 
tribution within the qubit enables the cancellation of (certain) 
decoherence modes. 

Contrary to spin-based quantum dot qubits, where 
decoherence-free subspaces can be created by combining 
quantum dots into logical units, charge-based qubits are much 
more difficult to isolate from the environment. In order to 
have decoherence-free subspaces for charge qubits one would 
need to restrict the operation to a subspace where charge is 
homogeneously distributed in space, no matter which basis 
states are chosen. However, this contradicts the very nature 
of a charge qubit (where readout depends on charge imbal- 
ance) and thus cannot be achieved. In our example of the 
three-dot qubit these facts become evident in the existence of 
two phonon modes that cannot cancel due to geometric con- 
straints inherent to the qubit. 

Decoherence can be mitigated in a number of other ways. 
For instance, for the three-dot qubit case we have studied, a 
substantial improvement with respect to the double-dot qubit 
can be achieved due to the lower frequency of operation and 
to an enhancement of the relaxation time by a factor of three. 

Another effective way to reduce the coupling to gapless 
bosonic modes is to choose a computational basis with a mul- 
tipole charge configuration. As we have shown, the multipole 
geometry attenuates the coupling to long wavelength acous- 
tic phonons by a factor proportional to a power law of the 
operation frequency. This power law grows rapidly with the 
multipole order. Thus, multipole configurations of charge can 
lead to quality factors enhanced by orders of magnitude in 
comparison to those obtained for double-dot qubits. How- 
ever, the effect is reversed at high frequencies of operation. 



The crossover frequency separating the two regimes is given 
by the inverse traversal time for a phonon to propagate across 
the qubit. For typical GaAs setups, this time is of the order 
of 30 ps (for dots 120 nm apart), indicating a crossover fre- 
quency in the range of 30 GHz. Since tunnel amplitudes usu- 
ally vary from tens to a few hundreds of /ieV, yielding quan- 
tum oscillations of about 2-20 GHz, there is a real advantage 
in moving toward multiple-dot qubits for current setups. How- 
ever, since phonons are not the only leading mechanism for 
decoherence in charge qubits, 15 as operation frequencies go 
down other sources of decoherence, not necessarily modeled 
by bosonic environments, may become dominant. In that case 
multiple-dot qubits might become less appealing. Investiga- 
tion of other decoherence mechanisms in multiple-dot qubits, 
such as electromagnetic environmental fluctuations and tele- 
graph noise due to charges trapped in the substrate, are the 
subject of ongoing investigations. 

Finally, we mention that a recent work has shown that 
pulse optimization is also a very effective way of minimizing 
the coupling to bosonic environments in quantum dot charge 
qubits.— 
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APPENDIX A: THE TWO-QUBIT REDUCED 
HAMILTONIAN 

The Hamiltonian of two three-dot qubits coupled by their 
bases [see Fig. |3b)] with inter-qubit couplings t' and t" has 
the following matrix form in the basis of Eqs. (fTTT>-(f20b (the 
lower off-diagonal block is omitted): 
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-t'e^/3 


-t"e-^/3 f'e 


-^/3-t"e^/3 


t"e^/3 
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Ui 
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Ui 
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V 











—v 
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Ui J 



(Al) 



[Note that the Hamiltonian for the tip-tip configuration is re- 
covered by setting t" = in Eq. ( IAlb .1 This Hamiltonian can 
be projected onto the two-qubit computational subspace by 
means of a Schrieffer- Wolff transformation. From Eq. (lAlb . 
we see that the Hamiltonian has the form ifjjff = Ho + Hi, 
where 



and T can be broken into two distinct 4x3 blocks, 



T = ( Ti T„ ) 



(A6) 



H = 




M 



Hi = 



T 
Tt 



(A2) 



and M and T are 6x6 and 4x6 matrices, respectively. Per- 
forming the Schrieffer- Wolff transformation and expanding to 
second order in H u 31 we get iJj^ff « iJ + (1/2) [5, Hi], 
where 



As a result, = -Ti B- 1 t} - T u B^ 1 T^. After some 

algebra, one finds that 



S = 





-l\r l T^ 



-TM- 1 




(A3) 



Thus the Hamiltonian has the block diagonal structure 



B' 




(A7) 



Vbasc 
Z I— II 



trrcd 

o 





M 



(A4) 



where Hfl d u 



= -TM- 1 Tl The matrix M can be broken 
into two identical 3x3 diagonal blocks, 



M = 



B 
B 



where u\ = (Ui + and U2 = (Ui — 2v) . The 

structure of B^ 1 can be substantially simplified by assum- 
ing v <§; Ui and neglecting v. This yields £f£fjj = 

(-1/C/i) (Tj T/ + T n T/j) . Carrying out the matrix multi- 



(A5) plications and setting t" = 0, we obtain Eq. ( 1231 . 
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